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Abstract. We present radial velocities for 29 B- and A-type stars in the field of the nearby association Perseus 
0B2. The velocities are derived from spectra obtained with AURELIE, via cross correlation with radial velocity 
standards matched as closely as possible in spectral type. The resulting accuracy is ~2-3 km s~^. We use these 
measurements, together with published values for a few other early-type stars, to study membership of the 
association. The mean radial velocity (and measured velocity dispersion) of Per OB2 is 23.5 ± 3.9 km s~^, and 
lies ~15 km away from the mean velocity of the local disk field stars. We identify a number of interlopers in 
the list of possible late-B- and A-type members which was based on Llipparcos parallaxes and proper motions, 
and discuss the colour-magnitude diagram of the association. 
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1. Introduction 

Perseus 0B2 is one of the nearest OB associations. It was 
discovered visually as a loose group of 15 bright, blue stars 
by Blaauw (1944), who subsequently confirmed his find- 
ing by means of radial velocity data from Moore's (1932) 
catalogue. Blaauw (1952; B52; see also Blaauw 1964) es- 
tablished membership for 17 bright O and B stars based 
on proper motion data. These include a number of well- 
studied spectroscopic binaries (e.g., AG Per), as well as 
the high-mass X-ray binary X Per. 

De Zeeuw et al. (1999; Z99) pubHshed an updated 
member list for Per 0B2, based on Hipparcos position, 
proper motion, and parallax data, containing 17 B- and 16 
A-type stars plus a small number of late-type stars. Based 
on extensive modeling of the kinematics of the Galactic 
disk, taking selection criteria in the Hipparcos Catalogue 
into account, Z99 concluded that a significant fraction of 
the A- and later-type stars identified as astrometric mem- 
bers are likely to be interlopers. This conclusion was con- 
firmed by Belikov et al. (2002a, b), who used Tycho-2 
proper motions as well as photometric information. These 
authors also identified nearly 1000 additional probable 
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members of Per 0B2 to 12th magnitude, and suggested 
that this association may consist of two subgroups. 

The Hipparcos satellite (ESA 1997) measured stellar 
positions, parallaxes, and proper motions, providing five- 
dimensional data in phase space. The sixth component, 
the radial velocity, is important, e.g., for improving mem- 
bership and expansion studies (e.g.. Brown, Dekker & de 
Zeeuw 1997). This is particularly true for Per 0B2, as 
the space motion of the association relative to the lo- 
cal disk population is mostly along the line-of-sight and 
proper motions of member stars are consequently small. 
Unfortunately, a homogeneous set of radial velocities is 
not available for the (early-type) Hipparcos members of 
Per 0B2. In this paper, we present new spectroscopic ob- 
servations of 29 B- and A-type stars in Per 0B2. 

While determining (relative) spectroscopic radial ve- 
locities for late-type stars using cross correlation tech- 
niques is feasible with precisions of tens of m s~^, sub- 
stantial effort is required to achieve even km s^^-level 
precisions for early-type stars (e.g., Verschuercn & David 
1999; Verschueren, David & Griffin 1999; Griffin, David & 
Verschueren 2000). Early- type optical spectra show few 
absorption lines and these lines are intrinsically broad 
(up to a few hundred km s^^). They are often broad- 
ened even further by stellar rotation and sometimes also 
show variability due to pulsations and/or stellar winds. 
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Fig. 1. Observing strategy. Each object exposure (either target or standard star) was preceded by a calibration se- 
quence, consisting of five bias exposures (Bl through B5), five flat field exposures (Fl through F5), and two wavelength 
calibration arc spectra (Wl and W2), and followed by an identical calibration sequence (B6 through BIO, F6 through 
FIO, W3 and W4). In order to minimize overhead, we used the trailing calibration sequence of each star also as the 
leading sequence for the next star. Time progresses from left to right. 



The resulting correlation peaks are broad, and can have 
significant substructure caused by the mixing of spectral 
lines (cf. figure 1 in Verschueren & David 1999). These 
effects complicate accurate centering of correlation peaks 
and thus the precise determination of (even relative) ra- 
dial velocities. The derived radial velocities of OBA-type 
stars, moreover, are known to depend on the spectral re- 
gion used in the cross correlation (e.g., Verschueren et al. 
1999). Furthermore, these stars emit their radiation pri- 
marily in the blue part of the spectrum, where CCD de- 
vices have reduced quantum efficiencies and slit-centering 
and atmospheric refraction errors are potentially signifi- 
cant (e.g., Verschueren et al. 1997). Finally, unrecognized 
multiplicity can be a significant source of error. 

We show here that, despite these significant com- 
plications, high-resolution spectroscopy combined with a 
careful observing strategy which includes measurement 
of many standard stars, does allow measurement of ra- 
dial velocities for B and A-type stars with an accuracy 
of a few km s~^. This makes it possible to improve sub- 
stantially the membership list for the early-type members 
of Per OB 2. We describe our observations in Sect. 2. In 
Sect. 3, we focus on the data reduction. The data analysis 
is presented in Sect. 4. The interpretation of the results 
and the conclusions are given in Sections 5 and 6, respec- 
tively. 

2. Observations 

The spectra presented here were obtained between Oct. 
31 and Nov. 13, 1997, with the 1.52 m Coude telescope 
at Observatoire de Haute-Provence (OHP). We used the 
AURELIE detector, a linear 2048-pixel photodiode array, 
which is relatively blue-sensitive (Gillet et al. 1994). We 
used grating R3 with a linear dispersion of ~16.5 A mm~^ 
and a resolving power R ^ 7000 in order to achieve a pre- 
cision of a few km in the final radial velocities; this 
precision is of the same order of magnitude as the typical 
error in tangential space motions due to Hipparcos proper 
motion errors at the distance of the Per 0B2 association 
(318 ± 27 pc; Z99), and comparable to the expected in- 
ternal velocity dispersion of the group. We used a Th-Ar 
lamp for wavelength calibration, as well as an image sheer 
with a circular 3" entrance aperture, which minimizes po- 



tential wavelength errors due to poor centering (Gillet et 
al. 1994). As our targets are hot B- and A-stars, which pri- 
marily show strong hydrogen and helium lines, we used the 
spectral range from 3800 to 4200 A. This choice resulted in 
the presence of about 10 higher-order Balmer and helium 
lines in a typical spectrum. The 3800-4200 A spectral re- 
gion also contains the Call H- and K-lines (K: A = 3933 A; 
H: A = 3968 A; cf. Figure EJ. In B-type stars, the K- 
line is mostly interstellar, but at spectral types around 
AO and later, this line, together with other metallic lines, 
becomes a notable feature intrinsic to the stellar spectra. 
For early-type stars, the Call H- and K-lincs thus trace 
gas along the line of sight instead of stellar kinematics 
(see, e.g., Sonncntrucker et al. 1999, who actually stud- 
ied the gas and dust distribution towards Per 0B2 using 
high- resolution spectra of bright stars in the association). 
Although the presence of (Doppler-shifted) H- and K-lines 
can potentially affect cross correlation results and the as- 
sociated radial velocities (Sect. 0J, the Call metal lines 
carry minimal weight in practice and their effects are neg- 
ligibly small. 

Z99 established membership of Per 0B2 by applying 
two independent selection methods (the convergent point 
method of de Bruijne [1999] and the Spaghetti method 
of Hoogerwerf & Aguilar [1999]) to all Hipparcos entries 
in a specific field on the sky (see table Al in Z99). This 
approach classified each star as either 'certain member' 
(acceptance by both methods), 'possible member' (accep- 
tance by one of the two methods exclusively), or 'non- 
member' (rejection by both methods). As targets for our 
observations, we selected the 33 B- and A-type certain 
members of Per 0B2 identified by Z99. We added a num- 
ber of their possible members, as well as some early- type 
members from B52 which were not confirmed by Z99. 

We observed the resulting list in order of decreasing 
brightness, ultimately obtaining high-quality spectra for 
29 distinct targets (Table^l. Exposure times ranged from 
10 to 30 min. Many targets were observed multiple times. 
In total, more than 7 nights of our 2-week observing cam- 
paign were weathered out completely. 

In addition, we repeatedly obtained high-quality spec- 
tra for 20 (candidate) radial- velocity standard stars. As 
an lAU-approved list of early-type radial velocity stan- 
dards does not exist, we selected these stars from various 
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Fig. 2. The stability of the detector, AURELIE. Each of the 131 star exposures has four associated wavelength 
calibration lamp spectra (Wl through W4): Wl and W2 were obtained before and W3 and W4 after each stellar 
exposure (Figure^. Lamp spectra W3 and W4 were generally taken ~30-45 min later than Wl and W2. Top panel: 
the filled circles show the shift (in km s^^) between the velocity scales inferred from Wl and from W2. This shift is, for 
our purposes, acceptably small (cr = 0.6 km s^^), with the exception of 21 cases (crossed symbols). These exposures 
are treated with extra care in this study. The dashed vertical lines separate different observing nights, the start dates 
of which are printed in the panel (year: 1997). Bottom panel: open circles and crosses denote the shifts between Wl 
and W3 and Wl and W4, respectively. Random shifts up to ~10 km s~^ occur (ct = 2.9 km s~^), most likely due to 
instability of the detector (Gillet et al. 1994). 



sources (Table El Sect. I4.l|l . In this selection, we tried to 
cover spectral type and luminosity class ranges as large as 
possible. We also preferentially selected stars with small 
rotation velocities. Although time consuming, the neces- 
sity to build up such a private standard star library, ob- 
served with the same instrumental setup as used for the 
target stars, has the advantage of working with a homo- 
geneous data set, and optimizes the accuracy of the final 
results (we use 'precision' for random and 'accuracy' for 
systematic errors). 



We obtained a total of 131 object exposures (target 
and standard stars together). We decided to take a cal- 
ibration sequence, consisting of five bias exposures, five 
flat fields, and two wavelength calibrations, before and af- 
ter every star was observed (Figure [T]l. In order to reduce 
overhead, we re-used the trailing calibration sequence of 
each star as the leading sequence for the next star to be 
observed. Although this procedure still resulted in a sig- 
nificant overhead, it helped to minimize the effect of in- 
strumental errors in our final results. 



3. Data reduction 

The data reduction was performed using IRAF software 
(Tody 1993). For each target exposure, we first subtracted 
the average bias exposure calculated from bias exposures 
Bl through BIO (Figure^]). In the next step, we similarly 
divided by an average flat fleld exposure. However, we 
rejected the first of the five flat field exposures in each 
calibration sequence (Fl and F6; Figure in order to 
allow for warming-up effects of the flat fleld lamp. We 
then manually removed cosmic rays from our data, after 
which we wavelength-calibrated the Th-Ar lamp spectra. 

In order to verify the stability of the detector, we 
cross correlated the four lamp spectra Wl through W4 
taken before and after each star was observed. We found 
generally good agreement between the wavelength scales 
inferred from the two lamp spectra taken either before 
(Wl and W2) or after (W3 and W4) the target expo- 
sure (cr = 0.6 and 0.7 km s^^, respectively; Figuresn]and 
121). However, cross correlating Wl with either W3 or W4, 
which were taken ~30-45 min later, showed random shifts 
of up to ±10 km s~^ (cr ~ 2.9 km s~^). These shifts effec- 
tively degrade the resolution of the detector, and are most 
likely due to thermal instabilities of the detector, which 
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Table 1. Observed target stars: B- and A-type members of Per 0B2 identified by Z99, supplemented with a number 
of early- type stars identified as members in pre-Hipparcos investigations (e.g., B52). The stars were observed in 
order of decreasing brightness, resulting in spectra for 29 distinct targets (62 exposures). Columns: (1) HD number 
(Hipparcos Catalogue field H71); (2) Hipparcos identifier (HI); (3) name; (4) V magnitude (H5); (5) spectral type 
and luminosity class; (6) number of spectra obtained (iV; subtractions refer to suspect exposures related to detector 
instability); (7) multiplicity (SB: spectroscopic binary; C: Hipparcos component binary — p and AHp in column 11; 
G: Hipparcos acceleration binary; S: Hipparcos suspected non-single — field H61); (8) Hipparcos membership (C: Z99 
certain member; P: Z99 possible member; B: Z99 non-member, but B52 member); (9) Hipparcos parallax tt (Hll; 
mas); (10) Hipparcos parallax error (Ttt (H16; mas); (11) Hipparcos components. 
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are known to be present at the level of ^ 0.1 pixel (hr ^) 
or -15 km s'^ (hr'^) (Gillet et al. 1994). Shifts be- 
tween W1/W2 and W3/W4 cannot be due to telescope 
re-pointings as AURELIE is located on a dedicated opti- 
cal bench in the telescope control room. We did not detect 
any systematic, long-term trend (spanning hours to days) 
in the detector zero point, consistent with the small re- 
peatability errors derived for a number of late-type stars 
that were observed throughout our campaign (Sect. 

Because the lamp spectra in the first calibration se- 
quence (Wl and W2) were taken immediately before the 
stellar exposure (Figure P), we decided not to use the 
two lamp spectra in the following sequence (W3 and 
W4), as these were taken after another five bias expo- 
sures and flat-field exposures. Unfortunately, even Wl and 
W2, taken within minutes of each other, show significant 



shifts (> cr = 0.6 km s^^) in 21 cases (crossed symbols 
in Figure 121). We reject the associated exposures in case 
of a standard star pointing (10 observations; see column 5 
in Table 121 and treat these exposures with special care in 
case of a target pointing (11 observations; see column 6 in 
Table CJ. 

In the final step of the data reduction procedure, 
we thus wavelength calibrated our spectra using the 
Wl Th-Ar lamp spectra. Examples of two fully reduced 
(continuum-subtracted) spectra are shown in Figure |2| 

4. Data analysis 

Radial velocities are generally determined by means of 
spectral cross correlation between the target and a suitable 
standard star (but see, e.g., Dravins, Lindegren & Madsen 
1999). Standard star spectra can either be obtained ob- 
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Fig. 3. Example of reduced (continuum-subtracted) spectra: HD 1438 (B8V; F = 6.10 mag; T — 30 min; left panel) 
and HD 24583 (B7V; V — 9.00 mag; T = 30 min; right panel). Note the near-absence of metal lines, the broad nature 
of the higher-order Balmer hydrogen lines and the helium lines, and the spectral differences in the helium lines. 



servationally or generated by means of stellar evolution- 
ary codes (e.g., Morse, Mathieu & Levine 1991). Whereas 
the latter approach has the advantage that the spectral 
type, luminosity class, rotation velocity, etc. of the tem- 
plate can be chosen representative of the target spectrum, 
significant uncertainties exist in the (atmospheric) models 
of early- type stars (e.g., Lennon, Dufton & Fitzsimmons 
1992). We chose to use our own standard star spectra for 
radial velocity calibration. Possible drawbacks in this ap- 
proach are, e.g., unrecognized multiplicity and errors in 
catalogued radial velocities and/or spectral classifications. 

4.1. Standard stars 

Our 20 early-type standard stars have been studied care- 
fully in the past, mostly with the aim of establishing bi- 
narity, and they have never shown signs of velocity vari- 
ability, although HD 27638, HD 196724, and HD 214994 
may have astrometric companions (Table|2J. Selection cri- 
teria in the construction of this standard star sample in- 
cluded visibility, visual magnitude, spectral type, luminos- 
ity class, and small rotational velocity. 

The published radial velocities for our standard stars 
generally agree to within ^1 km s^^, but differences up to 
~3 km s~^ occur. Some of these discrepancies may be due 
to differences in the absolute zero points used in differ- 
ent studies. In order to minimize the effect of this uncer- 
tainty, we preferentially selected the radial velocities for 
our standard stars from the source containing most mea- 
surements in absolute numbers (Morse et al. 1991; results 
based on CCD data). For the remaining stars, we used, in 
order of decreasing preference, Fekel (1985), Wolff (1978), 
Abt & Levy (1978), Gies & Bolton (1986), and Latham 
& Stefanik (1982). Exceptions to this rule were made for 
HD 196724 and 219188, for which we assume a radial ve- 
locity of —18.4 km s^^ (Fekel 1990; Liu, Janes & Bania 
1989) and +A8.0 km 8"^ (SIMBAD), respectively. Upon 
comparing the accordingly selected radial velocities with 
the values contained in the compilation catalogue of mean 
radial velocities of Barbier-Brossat & Figon (2000), we 



find a mean difference of 0.1 km s~^ and a standard devi- 
ation of 2.0 km s^^. These values (probably) mostly reflect 
zero-point differences (and random errors) and imply the 
accuracy and precision of our final radial velocities are 
^0.1 km s^^ and Jj2.0 km s~^, respectively. 

4.2. Cross correlation 

We cross correlated the spectra by means of the stan- 
dard IRAF task fxcor. In order to suppress noise, we 
applied high-frequency Fourier filtering to all spectra (cf. 
Verschueren & David 1999). We did not apply low-pass fil- 
tering, and used a Gaussian fitting function for centering 
of the correlation peak. This function was found empiri- 
cally to be a suitable choice. The detector instability did 
not justify a more refined fitting function. In order to sup- 
press the effect of continuum mismatch, we subtracted the 
continuum from both the template and the target spec- 
trum — and thus normalized both spectra to a continuum 
level of 1 — before each cross correlation. 



4.3. Repeatability errors 

We observed many of the objects, either target or standard 
stars, multiple times. Under the assumption that these 
stars are single, cross correlation of their spectra (on an 
object-by-object and exposure-by-exposure basis) should 
ideally result in vanishingly small shifts. Conversely, shifts 
measured in such an exercise provide a direct estimate of 
the random errors that will be present in the final radial 
velocities (e.g., Morse et al. 1991; Verschueren et al. 1997). 
In our case, these random errors include the temporal in- 
stability of the detector (Sect.lSJ. 

Figure 01 shows the repeatability errors inferred from 
our data, after excluding known spectroscopic binaries. 
We define the repeatability error for a pair of exposures 
of a given object as the relative shift (in km s~^) emerg- 
ing during their cross correlation. Sixteen of our standard 
stars and 18 of our Per 0B2 targets were observed multiple 
times, resulting in 86 repeatability errors (the two expo- 
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Table 2. Stars used as radial velocity standards. Columns: (1) HD number (Hipparcos Catalogue field H71); (2) 
Hipparcos identifier (HI); (3) V magnitude (H5); (4) spectral type and luminosity class (H76); (5) number of spectra 
obtained (iV; subtractions refer to rejected exposures due to detector instability); (6) multiplicity (G: Hipparcos 
acceleration binary; S: Hipparcos suspected non-single — field H61) (7) assumed radial velocity (km s^^; see Sect. 14. ll 
for details); (8) idem from Morse et al. (1991; table 6; table 4 for HD 23408); (9) radial velocity in km s^^ from Fekel 
(1985; table 3); (10) idem from Wolff (1978; range < 3 • standard deviation; tables 1 and 3); (11) idem from Abt & 
Levy (1978; table 1); (12) idem from Gies & Bolton (1986; table 3); (13) idem from Latham & Stefanik (1982; table 

1). 
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sures of the faint target HD 278942 lead to a repeatability 
error of ~16 km s~^; as these spectra are noise-dominated, 
we reject this object from this discussion). From Figure 21 
we conclude that the repeatability errors are at the level 
of ~0-3 km s~^, consistent with the expected noise due 
to thermal instabilities of the detector (Figure The 
mean repeatability error is o7 = 2.1 km s^"'^ (the median 
is 2.2 km s~^); this value could be slightly overestimated as 
a result of unrecognized multiplicity. A value of 2.1 km 
is comparable to repeatability errors quoted in the litera- 
ture (e.g., Morse et al. 1991: 1-3 km s~^; Verschueren et 
aL 1997: 0.7-1.4 km s'^). 



We also repeatedly observed two late-type stars 
(HD 18449, K2III, iV = 7; HD 219615, G7III, iV = 4). 
The repeatability errors for these stars are small, 0.6 and 
0.9 km s"""^, respectively, partly reflecting the relative ease 
with which late-type spectra can be cross correlated. The 
N — 7 exposures of HD 18449 were observed on 7 dif- 
ferent nights; the low value of the corresponding repeata- 
bility error is consistent with the absence of significant 
night-to-night instrumental zero-point shifts (Sect. |3| and 
Figure 13). 



4.4. Template (mis)match: standard star selection 

A radial velocity determination by means of cross corre- 
lation of an observed spectrum with that of a velocity 
standard can be biased if the spectrum of the target does 
not have exactly the same characteristics as that of the 
standard star. This template mismatch can, e.g., be due 
to differing rotational velocities (Sect. It can also re- 
sult from spectral type and/or luminosity class differences, 
which may translate into differences in atmospheric veloc- 
ity fields (convection, stellar wind, pulsation, . . .), in line 
blending, line blanketing, line asymmetries (due to Stark 
broadening), or in stellar continuum. 

We decided to use Hipparcos spectral classifications 
(catalogue fields H76 and H77) throughout this study. 
These values are extracted from ground-based compila- 
tions, mainly SIMBAD and the Michigan Spectral Survey 
(see ESA 1997 for details), and do not necessarily have 
a high and/or homogeneous quality. Classification errors 
(e.g., at subclass level) can thus not be excluded. In case 
the Hipparcos spectral classification lacked a luminos- 
ity class, we assumed the star was a dwarf. We ignored 
classification details/peculiarities such as 'p', 'n', and 'e'. 
Comparing the Hipparcos spectral types and luminosity 
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Fig. 5. Luminosity classes (LCI) and spectral types (SpT) of target stars (open circles) and standard stars (crosses). 
This diagram is used for selection of suitable standard stars for cross correlation. For a given target star with luminosity 
class LCI* and spectral type SpT,, we consider all Ns standard stars for which d < cJ^ax (eq.^. For spectral types later 
than B7 (>B7; dashed vertical line), we use dmax = 1-0 while for stars with SpT, = B7 or earlier, we use dmax = 1-5 
because of sparsity of (observed) standard stars. The circle demonstrates the selection of standard stars for a target 
star with SpT, = B8 and LCI* = V; for this object, we consider Ns — 3 velocity standards as appropriate templates 
(the three crosses falling within the circle of radius 1). The final radial velocity Uiad for this target star (eq.|31 Table |2Jl 
is obtained by means of a distance-weighted average of the individual radial velocities i^rad.i obtained using all Ng 
exposures of the i — 1, . . . , = 3 standard stars (Sect. 




5 
Repeatability error [km s"'] 

Fig. 4. Repeatability errors for 86 pairs of spectra (either 
target or standard stars; see Sect. 14.31 for details). In the 
remainder of this paper, we assume that repeatability er- 
rors (<Tr) are positive. The mean and median repeatability 
errors are 2.1 and 2.2 km s~^, respectively. 

classes with recent SIMBAD values and/or other litera- 
ture values did not reveal significant differences, except 
for the following target stars: HD 20987 (Hipparcos: B9 
^ Abt 1985: B2V), HD 23802 (Hipparcos: B9 => Guetter 
1977: B5Vn), HD 24012 (Hipparcos: B9 ^ Blaauw & 
van Albada 1963: B5), HD 24583 (Hipparcos: AO ^ 
Guetter 1977: B7V), and HD 278942 (Hipparcos: F2 ^ 
Cernis 1993: B3III). We discuss these stars individually in 
Sect. 14.81 We decided, for these five 'problem cases', to use 



the spectral classifications derived in the dedicated studies 
listed above. Alternatively, we could have used, e.g., the 
strength of the helium lines in our spectra as a coarse spec- 
tral type indicator (for example, the ratio Hel 4471/Mgn 
4481, not accessible from our data, is often used for de- 
termining the spectral types of stars; Jaschek & Jaschek 
1990). In case of line anomalies, line indicators may pro- 
vide misleading classifications. B-type stars, and notably 
the helium lines in their spectra, are known to be particu- 
larly sensitive to these problems. Dozens of peculiar stars, 
for example in the form of helium-strong and helium- weak 
stars, have been discovered in the nearby OB associations 
(Jaschek & Jaschek 1990), and we therefore decided not 
to use helium line strengths as spectral type indicators. 

In order to investigate the effect of spectral mismatch 
(including luminosity class mismatch), we performed a 
cross correlation of all possible combinations of all expo- 
sures of all 20 standard stars. After averaging over expo- 
sure pairs, we constructed a 20 x 20 'matrix' of average 
shifts between all pairs of standard stars. The diagonal ele- 
ments in this matrix, which is skew-symmetric, correspond 
to the repeatability errors fSect. After ordering the 
stars on spectral type, one can see that shifts grow when 
moving away from the diagonal. However, close to the di- 
agonal (i.e., within a few matrix elements, i.e., generally 
within a few sub-classes) , shifts are mostly relatively small 
(they are even generally insignificant in this exercise as a 
result of random errors). The results of this analysis indi- 
cate, nonetheless, that spectral mismatch can easily give 
rise to 5-10 km s~^ systematic errors or larger. For this 
reason, we decided, for a given target star exposure, not to 
use a single radial velocity standard as template. Instead, 
we used a set of standard stars with spectral types and lu- 
minosity classes 'similar to' the target, and averaged the 
corresponding radial velocities to one, final value. This has 
the additional advantage that radial velocity zero-point 
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Table 3. Radial velocities for 29 target stars. Each star occupies N lines, with N the number of spectra/exposures 
obtained (Table^. Columns: (1) HD number; (2) spectral type and luminosity class; (3) multiplicity (SB: spectroscopic 
binary; C: Hipparcos component binary; G: Hipparcos acceleration binary; S: Hipparcos suspected non-single — field 
H61); (4) the number of suitable standard stars Ng and, between brackets, the associated total number of exposures 
Nc (Sect. lT^ : (5) the average distance to the Ng standard stars (Sect. Eii| : (6) exposure number (1, ... ,7V; asterisks 
denote suspect exposures related to detector instability — see Sect. |31 for details); (7) heliocentric Julian date of the 
mid-point of the exposure (HJD 24507XX.XXXXX); (8) final distance-weighted radial velocity Wrad (km s~^; eq. |21 
in Sect. I4.4|l : (9) corresponding standard deviation iTt,,rad (km s~^; Sect. 14.6)) : (10) time-averaged radial velocity Wrad 
using the N exposures (km s^^; spectroscopic binaries show a dash; Sect. I4.6jl : (11) corresponding standard deviation 
o'v,rad (km s""'^; values smaller than 2.0 km s^^ are optimistic; Sect. I4.6|l : (12) literature radial velocity (km s~^); (13) 
source of literature radial velocity (column 12) and remarks (B52: Blaauw 1952; BvA: Blaauw & van Albada 1963; 
Z83: Zentelis 1983). 
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differences/errors in our standard star sample are washed 
out to some degree. 

In order to quantify the meaning of 'similar spectral 
types and luminosity classes', we defined a 'distance' d 
between any standard and any target star: 



(SpT, - SpT)2 + - . (LCU - LC1)2. 



(1) 



We then constructed a luminosity class-spectral type 
(LCl-SpT) diagram containing all target and standard 
stars (Figure 13 . From this diagram, we decided to define 
all standard stars for which d < dmax as suitable templates 
for a target star with luminosity class LCI, and spectral 
type SpT,. The unit of the distance d was chosen such 



that a difference of one spectral type subclass and a dif- 
ference of two luminosity classes both denote a distance 
of 1 (e.g., 'B8 - B6 -> d = 2' and 'V - III ^ d = 1'). The 
maximum distance (dmax) was chosen as small as possible, 
but still large enough for most target stars to have more 
than one associated standard star: 



1.5 
1.0 



for SpT, 
for SpT, 



< B7, 
> B7, 



(2) 



where <B7 means B7 or earlier and >B7 means later than 
B7. The radial velocity Vrad for the target is obtained by 
means of a distance-weighted average of the individual 
radial velocities Wrad,i obtained by cross correlating with 
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Table 3. Continued. 
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4.5. Template (mis)match: stellar rotation 



We carried out a simple experiment to investigate the ef- 
fect of rotational mismatch on the final radial velocities. 
We selected the US absorption line in the spectrum of the 
slow rotator HD 38899 (B9IV; vshii = 15-30 km s'^). 
We broadened this line, so as to mimic an increasing rota- 
tion velocity, with a Gaussian with increasing full-width 
at half maximum. We then cross correlated the broadened 
line with the original line and found no significant veloc- 
ity shifts (compared to the expected repeatability error of 
0^ — 2.1 km s^^), even for a broadening parameter which 
corresponded to a rotation velocity exceeding the break-up 
velocity. We repeated this exercise for HD 35708 (B2.5IV; 
vsmi — 10-24 km s^^) with similar results. Although 
this experiment is simplistic (it neglects, e.g., noise and 
effects due to changes in the limb darkening) , we conclude 
that rotational mismatch is not relevant at the level of 1- 
2 km s~^ (see also Verschueren & David 1999). Hence, our 



all No exposures of the i = 1, . . . , iVg suitable standard 
stars: 

where 

Wt = dmax - 0.5 ■ di. (4) 

This choice for Wi has the property that Wi{di = 
0)/wi{di = dmax) = 2, independent of spectral type. 

As an example. Figure |2| shows the spectra of HD 1438 
and HD 24583; the first star (B8V) has been used as a 
standard in the radial velocity determination of the latter 
(B7V). Although the signal-to-noise ratios and rotational 
velocities of the two spectra/stars differ, minor spectral 
differences are clearly visible, e.g., in the helium line at 
4026 A. The 'spectral-type-distance' d between the stars 
equals 1 (eq. ^ , which is the maximum allowed distance 
(dmax) for targets with spectral types B7V and later. 
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'bias' in preferentially selecting slowly rotating standards 
(Sect.|2) is not significant. 

4.6. Results 

Table 13 lists the results of applying the procedure de- 
scribed in Sect. 14.41 to the available target spectra. TablelHl 
also lists Ns and iVc, the number of standard stars and the 
number of standard star exposures used in determining 
Vrad, respectively, and the mean radial velocity for each 
exposure (eq. ; the corresponding standard deviation is 
denoted CTu^rad- The mean distance in luminosity class- 
spectral type space (Figure 13)) between the target and the 
Ne standard star exposures is denoted d = Nc~^ ■ X^i^i ^i- 
In general, the effect of template mismatch increases with 
increasing d. Table 13 also lists, for each target star, the 
time-averaged radial velocity TJrad, and the corresponding 
standard deviation ay^rad, of the N radial velocities corre- 
sponding to the N exposures obtained for each star. The 
quantity TJrad is the best estimate of the true radial veloc- 
ity of an object, provided it is single; Vrad is potentially 
less straightforward to interpret otherwise. 

Table 121 shows that, in general, <Jv,iad 3> cru.rad ^ ^7. 
The first inequality suggests that template mismatch, 
combined with possible errors in the assumed radial veloci- 
ties for the standard stars, is a significant error source, jus- 
tifying our averaging-approach. The similar magnitudes of 
<^v,Yad (with a mean value for all stars of 1.5 km s^^) and 
the mean repeatability error (a7 — 2.1 km s~^) provide a 
further indication of the validity of our approach. 

The errors tr^^rad on tJrad are sometimes remarkably 
small, down to a few tenths of a km s~^. Given that tem- 
plate mismatch can have an effect at the level of ~3.0~ 
5.0 km s^^ (corresponding to (Jv.rad), that the mean re- 
peatability error is 2.1 km s~^ (Sect. and that the 
expected error due to uncertainties in the radial velocities 
of the standard stars is also ~2.0 km s^^ fSect. E!T|l . we 
suspect that increasing cr^^rad to 2.0 km s^^ provides a 
more realistic estimate of the error. 

The 11 target exposures identified in Sect.Olas suspect 
as a result of detector instability are indicated with aster- 
isks in column 6 of Table |21 Comparing the associated ra- 
dial velocities with non-suspect exposures shows no signif- 
icant deviations, consistent with the relatively small mag- 
nitude of the effect established in Sect. 13 (0.6-2.9 km s~^; 
Figure 12). We therefore do not discriminate these suspect 
exposures further and treat them as normal in the remain- 
der of this manuscript. They thus also contribute to the 
time- averaged radial velocities I'rad- 

Table 0] shows, as a representative example, how the 
final, mean radial velocities Urad in Table |3I (eq. 13)1 are 
built up for the AO star HD 24970. This star was ob- 
served = 3 times and has Ng — 5 suitable standard 
stars (Figure O, giving rise to iVg = 5 radial velocities 
per exposure. The best-match radial velocity standard 
(HD 196724, AOV, d = 0.00, and w = 1.00) implies the 
target star's radial velocity is ~26 km s^^. The other four 



standards, however, which have rather closely matching 
spectral classifications, imply a velocity of '^20-25 km s""'^. 
Two explanations, which are not mutually exclusive, may 
be invoked to reconcile this discrepancy: either template 
mismatch explains the systematically lower values, or an 
error in the literature radial velocity of the AOV stan- 
dard HD 196724 explains the systematically higher value. 
Discriminating between these two options is impossible 
using the available data. The problem is, however, auto- 
matically 'solved' for/alleviated in the weighted-average 
approach adopted here. 

TableHshows, for the single AO stars HD 24970 {N = 
3; (Tr = 2.7 km s^^) and 23268 (A = 3; a, = 1.4 km s"^), 
for all combinations of the A^ = 3 object exposures and 
the As = 5 suitable standard stars, the standard devia- 
tions of the radial velocities corresponding to the different 
exposures of the standard stars (cf. TableEJ . The standard 
deviations are, in some cases, slightly larger than the ex- 
pected values (roughly the repeatability errors), suggest- 
ing that both errors in the literature radial velocities of 
standard stars and template mismatch contribute signifi- 
cantly to the total error budget. 

4. 7. Spectroscopic binaries 

Our sample contains five known spectroscopic binaries 
(SBs; a sixth SB, X Per/HD 24534, was classified as 
Per 0B2 member by B52, but was not observed by us; 
see Sect. l5.2?T|l . Three of these have previously been iden- 
tified as double-lined SBs (SB2s). The modest resolving 
power of the spectrograph-grating combination used by 
us {R 7000), combined with the very broad hydrogen 
absorption lines in the spectra of these SB2s, do not al- 
low a proper decomposition of the spectra in all of these 
cases (we did not optimize our instrumental setup for 
SB2s). We were therefore forced to treat these spectra 
as arising from a single star. The resulting radial ve- 
locities, derived by means of cross correlation using a 
template star matched to the spectral type and lumi- 
nosity class of the primary component, are of modest 
physical significance. Conceptually, they correspond to a 
luminosity-weighted average of the instantaneous radial 
velocities of the primary and secondary components, with 
template mismatch complicating matters (recall that the 
standard star was selected based on the primary compo- 
nent exclusively). While the composite spectrum we ob- 
served is a luminosity-weighted average of the instanta- 
neous Doppler-shifted spectra of the two components, the 
inferred radial velocity is not necessarily a luminosity- 
weighted average of the instantaneous radial velocities of 
the components, although a trend along these lines might 
be expected. 

The orbital elements from the literature for the five 
known SBs in our sample are listed in Table El We 
have one or a few measurements for each of these ob- 
jects. Repeat measurements were usually taken on the 
same night, typically within one hour, with the excep- 
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Table 4. Illustration of results underlying the mean radial velocities Vrad in Table |31(cf. eq. The star HD 24970 
(AO) was observed = 3 times and has — 5 suitable standard stars, giving rise to 5 radial velocities per exposure. 
These values are listed (in km s^^), together with the spectral classifications and distances d and weights w (see eqs^ 
and 0)) of the standard stars. The distance- weighted mean radial velocities with their errors are listed in the final 
column (in km s~^). Systematic differences between columns in this table can be due to template mismatch and/or 
incorrect literature radial velocities; in general, systematic differences between lines in this table can be due to random 
errors (statistics), instrumental shifts, and/or multiplicity (not likely in this specific case). 



HD 


10982 


196724 


196821 


27638 


58142 


"rad 


SpT 


B9V 


AOV 


AOIII 


B9V 


AlV 




d = 


1.00 


0.00 


1.00 


1.00 


1.00 




■w = 


0.50 


1.00 


0.50 


0.50 


0.50 




Exposure 1* 


26.2 


27.8 


24.3 


21.4 


22.8 


25.3 ± 3.3 


Exposure 2 


24.1 


26.1 


21.6 


20.1 


20.2 


23.2 ±3.5 


Exposure 3 


21.6 


23.8 


19.6 


16.9 


18.4 


21.0 ±3.5 



Table 5. The standard deviations of the radial velocities corresponding to the A^c different exposures of the Ng = 5 
suitable standard stars for all combinations of the different object exposures and suitable standard stars of the stars 
HD 24970 (AO; cf. TableU top part) and HD 23268 (AO; bottom part). The distance -weighted mean radial velocities 
with their errors are listed in the final column (in km s~^); time-averaged values, denoted Wrad, are provided in TableOl 
The quantity denotes the repeatability error (Sect. ^31). 



HD 


10982 


196724 


196821 


27638 


58142 




SpT 


B9V 


AOV 


AOIII 


B9V 


AlV 




d = 


1.00 


0.00 


1.00 


1.00 


1.00 




w = 


0.50 


1.00 


0.50 


0.50 


0.50 




N = 


2 


5 


3 


2 


4 




CTr = 


3.4 


2.3 


0.9 


0.2 


2.0 




Exposure 1* 


1.8 


3.4 


0.4 


0.2 


2.7 


±25.3 ±3.3 


Exposure 2 


2.3 


3.5 


0.9 


0.6 


2.7 


±23.2 ±3.5 


Exposure 3 


2.3 


3.5 


0.9 


0.7 


2.6 


±21.0 ±3.5 


Exposure 1 


2.1 


3.5 


0.8 


0.1 


2.6 


±3.8 ±3.6 


Exposure 2 


2.2 


3.5 


1.0 


0.6 


2.6 


±3.7 ±3.6 


Exposure 3 


2.3 


3.6 


1.0 


0.6 


2.6 


±4.8 ±3.6 



tion of o Per (Table EJ. As shown in Table El the ampli- 
tudes of the radial velocity variations are large, ranging 
up to ~180 km s~^ for AG Per. With the caveat men- 
tioned above for SB2s, our instantaneous measurements 
fall within these ranges in all five cases. 

For the two single-lined SBs (the SB Is), we used the 
known periods to investigate whether (i) our repeat mea- 
surements should have shown significant differences, and 
(ii) whether we could reconstruct what the expected Vrad 
should have been at the epoch of our observations. As the 
periods are known to sufficient accuracy, it is possible to 
do this, and we find agreement in both cases. We discuss 
each binary in more detail below. 

HD 23180 (o Per): This well-known SB2 (e.g., Stickland 
& Lloyd 1998) was observed twice (the first exposure 
being suspect), and then once more, ^^1.75 days later. 
According to Hipparcos, o Per is also an 1'.'019 visual bi- 
nary with AHp = 2.91 mag (Table ^ see Stickland & 
Lloyd for a discussion of this component). The Hipparcos 
data imply that our observations have primarily (~93% 
in flux) detected photons from the BlIII ± B2V SB. 



A simple-minded comparison of our radial velocities, 
obtained by simply ignoring the presence of the sec- 
ondary component, with predictions based on Stickland & 
Lloyd's (1998) orbital ephemeris shows the expected 'flux- 
weighted-average' trend (Table EJ. The third exposure of 
o Per shows some double helium lines with an average ve- 
locity separation of ~200 ± 30 km s~^; this is consistent 
with the expected velocity separation for this exposure of 
-230 km s-i (Table EJ. 

HD 23625: This object was studied by Blaauw & van 
Hoof (1963), who classified it as a B2V SB2. HD 23625 is 
also a Hipparcos component binary with AHp = 2.58 mag 
and p = 3'.'349 (Table^; these data imply that this fainter 
component cannot have been present in the 3"-diameter 
entrance pupil of the image slicer. We observed this object 
once, and find Wrad — 1-7 ± 3.0 km s^^ (our spectrum su- 
perficially looks single-lined, although a slight asymmetry 
in the line profiles is present; the radial velocity was again 
obtained by ignoring the presence of the secondary compo- 
nent in the spectrum) . According to Blaauw & van Hoof's 
ephemeris, we should have obtained —34.7 ± 19.1 km s~^ 
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Table 6. The five spectroscopic binaries (SBs) and their orbital elements. Columns for the first part of the table 
(literature data): (1) HD number; (2) periastron date To (HJD); (3) period P (days); (4) primary orbital semi-amplitude 
Ki (km s^^); (5) secondary orbital semi-amplitude K2 (km s~^; only for double-Hned spectroscopic binaries, i.e., 
SB2s); (6) systemic velocity 7 (km s~^); (7) orbital eccentricity e (see the footnote to the table for literature sources). 
Columns for the second part of the table (data from this study): (1) HD number; (2) exposure number (Table OJ; 
(3) instantaneous radial velocity (km s~^); (4) primary component prediction using the ephemeris provided above 
(km s^^); (5) secondary component prediction using the ephemeris provided above (km s~^; only for SB2s). 



HD 


To [HJD] 


P [days] 


Ki [km s^'] 


K2 [km s-'] 


7 [km s "'"] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


23180 


2445338.267 ±0.005 


4.4191447 ± 0.0000082 


111.8 ± 1.4 


155.0 ± 2.4 


+ 12.2 ±0.5 





23625 


2400000.0 ± 0.005 


1.940564 ±0.000004 


82.0 ± 2.0 


114.0 ±3.0 


±20.0 ±1.0 





24190 


2435401.00 ± 0.01 


26.111 ±0.001 


13.0 ±1.0 




+26.7 ± 1.0 


0.08 +0.01 


25799 


2435152.835 ±0.013 


0.9121679 ± 0.0000019 


17.6 ± 1.2 




+24.3 + 0.8 





25833 


2424946.5150 


2.0287298 


162.8 ± 1.1 


178.7 ± 1.2 


+24.7 + 0.9 


0.071 + 0.001 



HD 23180: Stickland & Lloyd (1998); 
HD 23625: Blaauw & van Hoof (1963); 
HD 24190: Lucy & Sweeney (1971); 

HD 25799: Morris et al. (1988); this is a non-radial pulsator; 
HD 25833: Gimenez & Clausen 1994 (cf. Popper 1974). 



HD 


Exposure 


Urad [km S ^] 


Wrad.l [km S ^] 


"rad,2 [km S ^1 


(1) 


(2) 


(3) 


(4) 


(5) 


23180 


1* 


-6.1 + 3.3 


-28.0+ 1.8 


67.9+ 2.5 




2 


6.4 + 3.1 


-21.6+ 1.8 


59.0+ 2.5 




3 


80.2 + 4.2 


108.7+ 1.6 


-121.6+ 2.5 


23625 


1 


1.7 + 3.0 


- 34.7+19.1 


96.1 + 26.5 


24190 


1 


21.3 + 4.6 


20.6+ 1.9 






2* 


22.2 + 4.6 


20.5+ 1.9 




25799 


1 


38.3 + 4.0 


40.3+ 2.1 






2 


38.4 + 3.9 


41.2+ 1.7 




25833 


1* 


31.6 + 1.2 


~190 


150 



for the primary component. This prediction agrees with 
our measurement at the 2tT-level. 

HD 24190: This object was identified as a low-amplitude 
B2V SBl by Blaauw & van Albada (1963; BvA). 
According to the improved ephemeris from Lucy & 
Sweeney (1971), our exposures (21.3 ± 4.6 and 22.2 ± 
4.6 km s~i) should have read 20.6 ± 1.9 and 20.5 ± 
1.9 km s~^, respectively. We conclude our measurements 
are in full agreement with Lucy & Sweeney's ephemeris. 

HD 25799: This object was classified as a low-amplitude 
eclipsing B3V SBl with a poorly determined/ variable pe- 
riod by BvA (P - 10.67 d and Ki = 20 km s'^). BvA's 
discovery of 'erratic velocity changes' was in fact the first 
indication that this object is not a normal binary. Morris 
et al. (1988) showed that an improved ephemeris (notably 
P — 0.9121679 d) could fit all previous spectroscopic data, 
but also showed, based on new photometric data, that the 
object is not a binary but a non-radial pulsator. This find- 
ing was confirmed later by Hipparcos photometry. We ob- 
served this object twice, obtaining v^ad = 38.3±4.0 km s~^ 
and 38.4 ± 3.9 km s~^. These observations are consistent 
with Morris et al.'s ephemeris within la, which predicts 
I'rad = 40.3 ± 2.1 and 41.2 ± 1.7 km s~^, respectively. 



HD 25833 (AG Per): This object is one of the few de- 
tached massive eclipsing SB2s, and is very well-studied 
(e.g.. Popper 1974; cf. Popper & HiU 1991). The sys- 
tem is eccentric (e ~ 0.071 ± 0.001) and shows apsidal 
motion with a period of 75.6 ± 0.6 yr (i.e., 4.76° yr~^). 
Spectroscopic ephemeris predictions are non trivial as a re- 
sult. As both components are roughly equally bright (V ^ 
7.46 and -7.88 mag; Table 6 in Gimenez & Clausen 1994) 
and have similar spectral types (B4V + B5; Gimenez 

6 Clausen), we expect the radial velocity inferred from 
our low-resolution composite spectrum to be near the 
systemic velocity, with a relatively small 'flux-weighted- 
average' oscillation due to orbital motion. Additional com- 
plications might arise from the presence of a third com- 
ponent, with Hipparcos parameters AHp = 1.81 mag 
and p = 0'.'803 (Table P). Gimenez & Clausen (1994) 
provide the most recent photometric eclipse ephemeris, 
from which we derive that the single, formally suspect, 
exposure we took of this object should have an associ- 
ated phase of —0.20 with respect to a secondary eclipse. 
This roughly implies radial velocities for the primary of 

7 + ~ +24.7 + 162.8 - +190 km s^^ and for the sec- 
ondary of 7-/^2 ^ +24.7-178.7-- -150 km s"\ so that 
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the expected velocity separation is ^340 km s"^. Not sur- 
prisingly, the hydrogen lines in our composite spectrum do 
not show clear double-lined signs (although an asymmetry 
of the line profiles seems present), but the much narrower 
helium lines are double. From these, we derive a mean ve- 
locity separation of 329 ± 17 km s~^ (based on 5 double 
lines), consistent with expectations. Ignoring the presence 
of the secondary component in the composite spectrum, 
we derive Wrad — 31.6 ± 1.2 km s^^ (based mainly on the 
broad hydrogen lines in the spectrum), which is indeed 
close to the systemic velocity. 

4.8. Comments on individual stars 

We noted in Sect. 14.41 that five targets have Hipparcos 
spectral types and luminosity classes which differ signif- 
icantly from other literature sources. In all these cases, 
we dropped the Hipparcos data and adopted the classi- 
fications from the dedicated studies. The implications of 
this choice for the inferred radial velocities are discussed 
below. 

HD 20987: The Hipparcos spectral type (B9) suggests 
Vra.d 40 km s-^ Abt (1985) and Roman (1978), in- 
dependently, quote a spectral classification of B2V, which 
suggests Wiad — —22.1 ± 0.8 km s^^ (Table |31 we used 

rfmax = 0.75). 

HD 23802: The Hipparcos spectral type (B9) suggests 
'^rad ^ —70 to —80 km s~^. Guetter (1977) quotes a 
spectral classification of B5Vn, which results in Vrad — 
-52.3 ± 6.2 km s'^ (Table El we used d^ax = 0.75). The 
three exposures give a relatively large spread in the in- 
ferred radial velocities (~12 km s~^), which, if significant, 
suggests that the radial velocity of this object is variable, 
most likely as a result of duplicity. 

HD 24012: The Hipparcos spectral type (B9) suggests 
Vrad ~ 19 km s^^. BvA quote a spectral classification of 
B5, which suggests frad = 26.8 ± 0.5 km s~^ (TableEJ. 

HD 24583: The Hipparcos spectral type (AO) suggests 
Viad 22 km s^^. Guetter (1977) quotes a spectral classi- 
fication of B7V, which suggests Vrad = 26.2 ± 5.7 km s~^ 
(Table EJ. Our five exposures seem to contain 1-2 out- 
liers (exposure 4 and, to a smaller extent, exposure 3). 
Excluding exposure 4 returns zJrad — 24.0 ± 3.0 km s~^; 
excluding both exposures returns Vi-ad ~ 25.4±1.6 km s^^. 

HD 278942: This faint object is a Hipparcos compo- 
nent binary {p = 0'.'149 and AHp = 0.78 mag; Tabled). 
An IRAS ring in the interstellar medium around this 
star explains the various colour and spectral type mea- 
surements reported in the literature (see Z99 for details; 
Hipparcos/Tycho: SpT = F2 and B - F = 1.13 mag; 
SIMBAD/AGK3: SpT = B5 and B - F = -0.1 mag; 
Cernis 1993: SpT = B3III; Andersson et al. 2000: SpT = 
O9.5V-B0V; the latter authors suggested that the IRAS 
ring, which is also visible at radio wavelengths, is an HII 
region associated with the object). Our spectra, although 
noisy, show that the object is a B-type rather than an F- 



type star. Wc therefore adopt the B3III spectral type from 
Cernis (1993), who identified HD 278942 as a 'possible 
photometric B3III + F5I binary'. This choice implies a ra- 
dial velocity of ^31 km s~^ (TableEl despite a repeatabil- 
ity error of ~16 km s~^, the inferred radial velocities of the 
two exposures are consistent at the level of '^l km s^^). 
Using the SIMBAD/AGK3 B5 spectral type would result 
in a radial velocity of ^^38 km s^^; Andersson's classifica- 
tion would imply a radial velocity of ^^70-90 km s^^. 

HD 281159: We observed this star twice (in subsequent 
exposures), obtaining v^ad — 9.1±1.8 and 7.8±2.1 km s~^. 
Both SIMBAD and Dufiot et al. (1995; the WEB cata- 
logue of radial velocities) list a literature radial velocity of 
-|-14 ± 5 km s^^, which seems in reasonable agreement at 
first sight. There is, however, a long history behind these 
values. B52, based on Moore's (1932) results, lists a radial 
velocity of +32 km s^^ and sets a radial velocity variabil- 
ity flag. Wilson & Joy (1952) list '-I-25 km s^^ (variable 
radial velocity)'. These authors based their mean value 
on 4 measurements: -1-42, -1-20, -1-80, and —43 km s~^. 
The General Catalogue of Radial Velocities (GCRV; 1953) 
lists +25 ± 5 km s^^ and does not mention velocity vari- 
ability. Petrie & Pearce (1961) list '+7.0 ± 4.7 km s"^ 
(variable radial velocity)'. These authors based their mean 
value on 6 measurements: +34, —10, —3, —4, +19, and 
+16 km s-i. Evans (1967; 'The revision of the GCRV) 
lists +14±5 km s^^ based on 10 measurements. This mean 
value is the number-of-measurements- weighted average of 
the mean results of Wilson & Joy (4 measurements) and 
Petrie & Pearce (6 measurements). Evans did not copy 
any of the variability flags in the literature, and his re- 
sult (+14 ±5 km s~^) erroneously suggests that the radial 
of HD 281159 is well-defined and reliable. Evans's results 
were copied by both Dufiot et al. (1995) and SIMBAD. 
We thus conclude that the agreement between the results 
of the latter two literature sources with our measurements 
is spurious and we suspect that this object is an SB. As 
the radial velocity determinations by Wilson & Joy and 
Petrie & Pearce have been published without the asso- 
ciated epochs, we cannot analyse the data for this star 
in the light of an SB model. Our mean radial velocity 
{vrad = 8.4 ± 0.9 km s~-^) should be treated with care. 

4.9. External accuracy check 

Our list of 29 stars with v^ad measurements contains the 
five confirmed spectroscopic binaries discussed in Sect. 14. 71 
and we have seen that our measured values are consis- 
tent with earlier measurements. We were able to find re- 
ported radial velocity measurements for 15 of the remain- 
ing 24 objects. The main sources, in order of decreas- 
ing preference, are Zentelis (1983; Z83), Blaauw & van 
Albada (1963; BvA), Blaauw (1952; B52), and Grenier et 
al. (1999). 

Figure shows our measured values (wqhp) for these 
stars versus those in the literature {vm ) . We have excluded 
HD 23268. Its literature radial velocity from Dufiot et al. 
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(Sect. I5.3II . and then analyse the colour- magnitude dia- 
gram (Sect.Eir 



20 

Radial velocity [OHP; km s"'] 

Fig. 6. Comparison of OHP and hterature radial veloci- 
ties. The filled circles denote 15 single stars (i.e., non-SBs). 
One of these, HD 23268 with OHP and literature radial 
velocities of 4.1 ± 0.6 and -20 ± 8.5 km s^^ (Duflot et al. 
1995), respectively, falls outside the plot. The four open 
squares refer to the 2-1-2 instantaneous measurements of 



the two single-hned SBs HD 24190 (wrad 
and 25799 (i^rad 38 km s"!). 



22 km s" 



(1995) of —20 ± 8.5 km s~^ is based on 3 measurements 
with unknown source; it is consistent with our measure- 
ment of 4.1 ± 0.6 km s~^ at the 3a level. There is rea- 
sonable agreement between the literature values and our 
measurements. A small offset from the diagonal can be ex- 
pected, as we have set the zero of our scale by the choice 
of standards, while in literature studies this is done in 
various different ways (cf. Sect. I4.1|l . We find that the 
weighted mean and dispersion of A = wqhp ~ vut are 
0.5 and 3.3 km s~^ (the straight mean and dispersion of 

A = (t;oHP-fiit)/-y/o'^,oHP +o'u,iit are 0.1 ±0.9, compared 
to the expected value of ± 1). These values are for the 
complete set of single stars plus the 2+2 SBl exposures for 
HD 24190 and 25799 (Figure EJl. We excluded HD 23268 
and the three 'outliers above the diagonal' (HD 23060, 
23478, and 24012) from this sample. Removing the SBl 
exposures from the sample results in 0.9 and 3.5 km s"^ 
(and 0.2 ± 1.0), while adding the three outliers results in 
-1.5 and 5.4 km s'^ (and -0.4 ± 1.6). 



5. Interpretation 

We now use our measurements, together with those for 
other certain and proposed members of Per 0B2, to study 
the association. We first discuss the distribution of ra- 
dial velocities (Sect. I5.1|l and membership of individual 
stars (Sect. I5.2|l . briefiy address the internal structure 



5.1. Distribution of radial velocities 

We use our own measurements for 24 stars, and use the 
reported systemic (7) velocities for the five spectroscopic 
binaries (Table in our sample of candidate Per 0B2 
members. The results are plotted in Figure [3 which shows 
the distribution of these objects on the sky, flanked by 
plots of Vrad vcrsus galactic longitude £ and latitude b. We 
use different symbols for our new measurements (squares) 
and the literature values (triangles). 

The Vrad versus i plot clearly shows a clump of stars 
near Vrad ^ 23 km s^^, containing many of the classical 
Per 0B2 members (B52). The tirad versus b plot shows 
a very similar separation. In addition to a few outliers, 
there is a second clump of stars, covering a larger range in 
£, with a small dispersion. The measured values coincide 
very nicely with the expected Wrad for field stars in the 
direction of Per 0B2, i.e., those obeying simple galactic 
rotation (in this direction, however, most of the observed 
radial velocity is reflected Solar motion) . We conclude that 
these objects are unrelated field stars. 

The radial velocity separation between the association 
and the Galactic disk allows us to determine the mean ra- 
dial velocity of the group. Our list of 29 stars contains 19 
stars with radial velocities between 10 and 35 km s~^. 
From this list, we reject the likely non-member o Per 
(HD 23180; see Sect.ESl- Figure El shows the radial veloc- 
ity histogram of the 18 (candidate) members. From this 
histogram, we obtain a mean radial velocity of 23.5 km s""'^ 
and an associated dispersion a = 3.9 km s^^. This mean 
velocity is consistent with the value derived by Blaauw 
(1944), 19.4 ± 1.7 km a"^. 

The dispersion of 3.9 km s^^ amongst the measured 
radial velocities in Per 0B2 provides an a posteriori ex- 
ternal check of the accuracy of our measurements, as it is 
an upper limit on this (the measured dispersion arises from 
measurement errors, internal velocity dispersion in the as- 
sociation, unrecognized duplicity, and/or the presence of 
non-members in the sample). The internal dispersion is 
probably only 1-3 km s^^ (Z99), which suggests that the 
external accuracy is ^3 km s^^. This is in harmony with 
the repeatability errors and the analysis of Sect. 14.91 



5.2. Membership 

The clean separation of the association and the local disk 
stars in radial velocity (Figured) makes it possible to im- 
prove some of the earlier membership assignments (no- 
tably B52 and Z99) . We discuss these in some detail here. 



5.2.1. Blaauw 1952 (B52) 

B52 identified 17 stars as member of Per 0B2. Z99 con- 
firmed eight of these: HIP 17313, 17735, 18081, 18111, 
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18246, 19039, 19178, and 19201 (HD 22951, 23625, 24131, 
24190, 24398, 25539, 25799, and 25833). Four of these 
eight are SBs (Table IHJl; their systemic radial velocities 
(from the literature) cluster around 20-25 km s~^. The 
OHP radial velocities of the other four objects cluster in 
the same range. We thus conclude that all eight B52/Z99 
Per 0B2 members are genuine members. 

Z99 identified four B52 members as possible astro- 
metric member: HIP 17387, 17448, 17465, and 18434 
(HD 23060, 23180, 281159, and 24640). We observed all of 
them, and conclude that HD 23180 (o Per), a SB with a 
systemic velocity of 12.2 km s~^ fSect. I^77|l . is most likely 
a non- member. According to their spectral types and ra- 
dial velocities, HD 23060 and 24640 are most likely mem- 



bers. HD 281159 is possibly a SB and could be a member 
(Sect. EH). 

The remaining five B52 members are — by definition 
— Z99 non-members: HIP 16203, 16518, 17631, 18350, 
and 18614 (HD 21483, 21856, 23478, 24534, and 24912). 

Both HD 21483 and 21856 were rejected by Z99. From 
the radial velocities, we conclude that HD 21483 is clearly 
a non-member (this object also has a deviating position on 
the sky). However, HD 21856 could be a member, certainly 
given its early BIV spectral type, although its parallax 
is 'small' (1.99 ± 0.82 mas), its radial velocity is 'large' 
(29.7 ± 3.3 km s~^), and its (Hipparcos and Tycho-2) 
proper motion deviates from the mean of Per 0B2. 
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Fig. 8. Radial velocity histogram of all 18 stars with OHP 
radial velocities between 10 and 35 km (we excluded 
o Per/HD 23180; Sect.ESi- We find a mean radial velocity 
of 23.5 km s"'^ (vertical line) and an associated dispersion 
of 3.9 km s~^ (dashed vertical lines). The curve denotes 
a reference Gaussian with a mean of 23.5 km s^^ and 
tr = 3.9 km s~^. See Sect. l5.ll for details. 

HD 23478 (B3IV...) was not tested for membership in 
Z99 because the Hipparcos astrometric data quality in- 
dicator was large (H30 = 3.21). Its OHP radial velocity 
(15.8 ± 4.0 km s'^), its parallax (4.19 ± 1.03 mas), and 
its Hipparcos and Tycho-2 proper motions are consistent 
with membership. 

HD 24912 is ^ Per, a celebrated run-away star (Blaauw 
1961; Hoogerwerf, de Bruijne & de Zeeuw 2001). It moves 
away from Per 0B2 with a relative radial velocity of 
~40 km s~^; we did not observe this object. 

The last of the remaining B52 members is X Per 
(HD 24534). This 09.5pe high-mass X-ray binary was 
identified as member by Blaauw (1944, 1952), although its 
absolute magnitude gave rise to doubts on this classifica- 
tion. Z99 did not consider it because the Hipparcos proper 
motion was of insufficient quality, and for this reason we 
did not observe it. The reported systemic radial velocity is 
uncertain at the level of 50 km s~^ (B52; Wackerling 1972; 
Hutchings et al. 1975; Duflot et al. 1995), so membership 
remains uncertain. 

5.2.2. De Zeeuw et al. 1999 (Z99) 

We observed 14 stars in addition to 15 of the 17 B52 
members (Sect. 15. 2. If) . All of these are certain or possi- 
ble astrometric Per 0B2 members from Z99 (see Sect.|2). 
Based on their OHP radial velocities, we conclude that 
HD 18830, 19216, 19567, 20113, 22114, and 23268 are un- 



related field stars in the disk ('interlopers'). Based on their 
OHP radial velocities and the Hipparcos data, we conclude 
that HD 24012, 24583, 24970, and 26499 are members of 
Per 0B2. We briefly discuss the remaining four stars be- 
low. 

HD 20987: The OHP radial velocity (-22.0 km s^^) sug- 
gests that this Hipparcos acceleration binary is a back- 
ground field star. This suggestion is strengthened by (i) 
the star's position on the sky '^10° away from the main 
body of the association; (ii) its relatively small parallax 
(1.80 ± 1.08 mas); (iii) its early spectral type (B2V) com- 
bined with its relatively faint magnitude {V = 7.87 mag). 

HD 278942: This faint object is a Hipparcos component 
binary (TableU)!. Cernis (1993) identified HD 278942 as a 
'possible photometric B3III -I- F5I binary'. Its presence in 
the Per 0B2 cloud and its parallax (4.83±1.21 mas) make 
it clear that the object is roughly at the same distance 
as Per 0B2 (-318 pc; Z99). We find a radial velocity of 
-38 km s^i for a B5 spectral type (SIMBAD/AGK3) and 
of —31 km s^^ for a B3III spectral type (Cernis 1993; see 
Sect. 14.81 for details). The uncertain spectral type, com- 
bined with duplicity, introduces a relatively large uncer- 
tainty in the inferred radial velocity. We suspect that this 
object belongs to Per 0B2. 

HD 23597: Although the OHP radial velocity is small 
(—16 km s~^), there is no compelling reason to believe it 
is not a member of Per 0B2. 

HD 23802: This B5Vn star is a Hipparcos component 
binary with a parallax of 3.09 ± 1.21 mas. At first sight, 
the OHP radial velocity (——50 km s^^) seems inconsis- 
tent with membership. The three OHP exposures, how- 
ever, give a relatively large spread in the inferred radial 
velocity (—12 km s~^), which, if significant, might be in- 
terpreted as resulting from duplicity. This in turn could 
imply that the object does belong to Per 0B2. Indirect 
evidence in favour of the above arguments is provided by 
the highly negative value of the OHP radial velocity: the 
correspondingly inferred distance of the object (>4 kpc) 
is inconsistent with its spectral type and magnitude. 

5.2.3. Literature radial velocities 

The analysis in Sections I5.2.1ll5~2.2l leaves three Z99 cer- 
tain members which we did not observe, but for which 
radial velocities are available from SIMBAD and BvA. 
The measurement for HD 23244 (AO) is highly uncertain 
(w,^d = 15 ± 12.4 km s-i from SIMBAD), but is for- 
mally consistent with membership. HD 20825 (G5III) is a 
clear non-member at 6.2 ± 1.2 km s'^ HD 281157 (B3V; 
a Hipparcos component binary) has a radial velocity of 
20.9 ± 2.1 km s^^ from BvA, and hence is a member. 

5.2.4. Conclusion on membership 

It follows from Sections 15. 2. II 15"!^^ that some of the clas- 
sical/possible members are in fact field stars. This shows 
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that the proper motion selection, while good, is not suffi- 
cient (this statement is particularly true for Perseus 0B2 
as the relative motion of the group is mostly along the 
radial direction/line of sight). Indeed, Z99 estimated the 
fraction of interlopers in their list of Hipparcos members, 
and concluded that for Per 0B2 at 318 pc, this is 70-100% 
over the range B9, A and later. These numbers agree with 
the Tycho-2 membership analysis of Belikov et al. (2002b) 
and also with the findings of this investigation. 

Table [7| summarizes our membership assignment for 
Per 0B2, based on proper motions (B52, Z99), radial ve- 
locities (mostly this work), and photometry (Sect. 15.4(1 . 
Membership derived from radial velocity data (col- 
umn 5) is based on a ±2tT criterion. The remaining 
stars are, by definition, classified as non-members, with 
the exception of the following five 'special cases', all of 
which are provisionally classified as possible members: 
HD 23802 and 281159, two suspected (spectroscopic) bi- 
naries, HD 278942, a faint target with an uncertain spec- 
tral type and low-quality spectra, and HD 24534 and 
23244, which have highly uncertain literature radial veloc- 
ities. Final membership (column 6) is based on combin- 
ing Hipparcos/astrometric (Z99) and radial velocity mem- 
bership (columns 4 and 5) following the logical scheme 
Col. A + Col. 5 — > Col. 6, where C + C — > C, P + P — > 
P, N + N — > N, C + P — > P, C + N — > N (as- 
trometric interloper), P + C — > C (astrometric binary), 
P + N — > N, N + C — > C (astrometric binary), 
N + P — > N (combination not present); C, N, and P 
stand for certain, non-, and possible member. In practice, 
final membership in our sample is effectively the same as 
radial velocity membership. 



5.3. Internal structure 

B52 did not distinguish subgroups for Per 0B2, although 
he did find subgroups for several other OB associations 
(cf. Blaauw 1964). However, Blaauw speculated that sub- 
groups would be created as the Per 0B2 association would 
evolve. Mirzoyan et al. (1999) claimed to have found sub- 
structure and expansion in Per 0B2 with the Hipparcos 
analysis of 17 bright stars in the association. Mirzoyan, 
using Hipparcos Input Catalogue (HIC) radial velocities, 
found two 'subgroups': one centered around -1-17.4 km s^^ 
and one around -1-26 km s~^. Belikov et al. (2002b) also 
presented evidence for two subgroups. We see no evidence 
for subdivision of Per 0B2 in Figure |S1 but our sample 
is small. For the particular case of Per 0B2, Hipparcos 
astrometric radial velocities offer no viable alternative to 
study internal structure issues: as the relative motion of 
the association is mostly along the line of sight, Hipparcos 
astrometric radial velocities are expected to be accurate to 
^^2.5 km s~^ (Dravins et al. 1999), similar to the accuracy 
of the spectroscopic radial velocities presented here. 



5.4. Colour-magnitude diagram 

Figure El shows the observed and de-reddened colour- 
magnitude diagrams (CMDs). In order to construct these, 
we removed the eight non-members identified above the 
horizontal line in Table [7| (denoted by 'N' in column 6), 
but retained HD 24970 as this is a possible member based 
on the astrometry and radial velocity measurement. We 
added the two classical members from B52 that we did 
not observe (^ Per/HD 24912, which is a 'member' in some 
sense, and X Per/HD 24534, which is a doubtful member; 
Sect. I5.2"T)| as well as HD 281157 which we classified as 
member based on a literature radial velocity fSect. [5?T^ . 
We retrieved Johnson B and V magnitudes for the result- 
ing 24 stars and derived individual extinctions Ay from 
both Stromgren and Johnson photometry using Ry = 3.2. 
We adopted Ay = 1 mag, roughly the mean extinction for 
Per 0B2 (e.g., Belikov et al. 2002b), for the star HD 26499 
which lacked appropriate photometry; this star is labeled 
'4' in Figure El 

The right panel of Figure El displays a tight main se- 
quence. The stars HD 278942 and 281157 (labels '1' and 
'2', respectively) both suffer from large extinctions {Ay — 
4.75 and 2.93 mag, respectively). The accuracy of our 
simple extinction correction is such that the de-reddened 
CMD locations are fully consistent with membership. The 
locations of two B5V stars with 'peculiar'/deviating ra- 
dial velocities (HD 281159 with Ujad — 8.4 km s~^ and 
HD 23802 with Vra.d = -52.3 km s-\ labels '3' and '6', 
respectively) are also consistent with membership, which 
suggests these stars could be multiple. Our adoption of the 
average extinction to Per 0B2 for HD 26499 (label '4', see 
above) puts it indeed on the dereddened main sequence. 

The CMD location of HD 24970 (label '5') raises doubt 
about its physical association with Per 0B2. Although the 
observed radial velocity (23.2 ± 2.2 km s^^) is consistent 
with membership, we suspect that, given the relatively 
bright magnitude for its spectral type (AO), its relatively 
'large' Hipparcos parallax (4.95 ± 0.99 mas), and given 
the fact that this star is not a certain but 'only' a pos- 
sible Z99 member, it is in fact a foreground object (we 
found Ay = 0.64 mag, roughly the foreground extinction 
to Per 0B2; Rydgren 1971). The peculiar location of X Per 
can be ascribed to the emission lines in its spectrum. The 
locations of the runaway star ^ Per (label '8') and the 
supergiant ( Per (label '10') are consistent with member- 
ship. Finally, o Per (label '9') seems to fit the zero-age 
main sequence, even though its radial velocity indicates 
it is a non-member (the object is, nonetheless, located in 
the Perseus molecular cloud). 

While our radial velocity study has removed some in- 
terlopers from earlier member lists of Per 0B2, including 
HD 21483 (B3III) and o Per (BlIII), most of these have 
mid-B spectral types or later. As a result, the location of 
the main sequence turn off is unchanged, and so are the 
age determinations based on this. 
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6. Concluding remarks 

We have presented radial velocity measurements derived 
from high-resolution spectroscopy for 29 candidate B- and 
A-type members of the association Per 0B2 obtained with 
AURELIE at Observatoire de Haute-Provence. The radial 
velocities were found via cross correlation with the spectra 
of observed standard stars with known radial velocities. 
A careful observing strategy and data reduction, rigor- 
ous removal of suspect exposures caused by an unstable 
detector, and use of multiple standard stars to minimize 
template mismatch, resulted in derived radial velocities 
accurate to ^^2-3 km s~^. We have checked our error es- 
timates by comparison with literature measurements. 

The new measurements, together with those for some 
additional candidate members from earlier studies, show 
that Per 0B2 lies offset by ^^15 km from the Wrad dis- 
tribution for field stars in this direction. The dispersion 
amongst the measurements for the Per 0B2 members is 
only 3.9 km s~^, which provides an upper limit to our 
radial velocity accuracy. The internal one-dimensional ve- 
locity dispersion of the association is smaller than we can 
measure, but is likely to be ~l-3 km s'-^. 



The radial velocity measurements confirm many of the 
classical astrometric members, but not all. Some earlier 
'secure' members turn out to have discrepant Vmd val- 
ues; this means they are either previously unknown spec- 
troscopic binaries or are genuine non-members. The bulk 
of these discrepant stars have Urad of about 6 km s~^, 
which puts them right at the velocity of the disk (field) 
stars in this direction, suggesting they are non-members. 
Removing these from the sample results in a very nar- 
row main sequence in the colour-magnitude diagram. The 
two outliers in this diagram are the X-ray binary X Per, 
whose membership remains uncertain, and the AO star HD 
24970. Even though its radial velocity and proper motion 
are consistent with the space motion of Per 0B2, we con- 
clude HD 24970 is a foreground object. 

The classical study by BvA had already shown 
that a significant fraction of the B-type members of 
Per 0B2 are spectroscopic binaries. Our study has identi- 
fied HD 23802, HD 281159, and possibly also the peculiar 
object HD 278942, as additional candidate spectroscopic 
binaries. It would be worthwhile to obtain multi-epoch 
spectroscopy for all the B-type members of PerOB2, in 
order to derive accurate orbits for the spectroscopic bi- 
naries, and possibly identify additional ones. The effort 
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required is significant, but would be invaluable in deter- 
mining the present day mass function of the association, 
and the nature of the binary population. 

The modest accuracy of the Hipparcos parallaxes 
and proper motions for groups such as Per 0B2, com- 
bined with the relatively bright limiting magnitude of the 
Hipparcos Catalogue, did not allow reliable identification 
of members later than ^B9 AO, even though these are no 
doubt present (see, e.g., Preibisch et al.'s 2002 study of 
the Sco-Cen OB association). A first step to probe mem- 
bership of Per 0B2 for later spectral types was made by 
Belikov et al. (2002a, b), who constructed a list of ~1000 
possible members of Per 0B2 based on Tycho-2 proper 
motions and photometric information. Our results demon- 
strate that a natural next step to distinguish members 
from interlopers in this list would be to obtain spectra 
from which to derive radial velocities (as well as spectral 
types and luminosity classes). 

ESA's future astrometry satellite GAIA (Perryman et 
al. 2001) will extend proper motion studies to much fainter 
limiting magnitude than Hipparcos, and will allow kine- 
matic identification of low-mass members, including all 
pre-main sequence objects, in Per 0B2 and other nearby 
associations. GAIA will also allow study of more distant 
groups in unprecedented detail. Besides astrometry, GAIA 
will also collect multi-epoch spectroscopic measurements 
for all stars down to F ~ 17 mag. However, due to the 
specifics of the wavelength range used (^^250 A around 
the Call-triplet near 860 nm), the resulting radial veloci- 
ties for early-type stars will likely be of modest precision. 
The currently planned ground-based all-sky multi-object 
spectroscopic radial-velocity survey RAVE uses the same 
spectral range, and may hence suffer from the same limita- 
tion. Our study shows that dedicated (multi-epoch) spec- 
troscopic studies of the early-type stars in the nearby asso- 
ciations are feasible, and are crucial in determining mem- 
bership and the nature of the population of high-mass 
stars in these young stellar groups. 
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Table 7. Membership and photometry of Per 0B2. We 
first list the 29 objects observed by us (of. Table 
Below the horizontal line, we added two classical mem- 
bers from B52 that were not observed (Sect. 15.2. l|l and 
three members from Z99 for which literature radial ve- 
locities exist (Sect. IF"2.3f) . Columns: (1) HD number; (2) 
spectral type; (3) classical membership assignment (B52; 
C: certain member; P: possible member; N: non-member); 
(4) Hipparcos membership (Z99); (5) membership based 
on radial velocity (mostly this work); (6) final member- 
ship, based on columns 4 and 5; in the case labeled with 
an asterisk, also photometric information was used; (7) 
Johnson B magnitude (from Tycho); (8) Johnson V mag- 
nitude (from Tycho); and (9) visual extinction Ay from 
Stromgren or Johnson photometry (Sect. I5.4|l . We do 
not list photometric data for eight clear non-members 
(Sect. 15. 4|) . Membership for the run-away star ^ Per 
(HD 24912) comes down to semantics; we have indicated 
it by the label R. 
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